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The small- a; deep inelastic scattering in the saturation region is governed by the non-linear evo- 
lution of Wilson-lines operators. In the leading logarithmic approximation it is given by the BK 
equation for the evolution of color dipoles. In the NLO the nonlinear equation gets contributions 
from quark and gluon loops. In this paper I calculate the quark-loop contribution to small-x evolu- 
tion of Wilson lines in the NLO. It turns out that there are no new operators at the one- loop level 
- just as at the tree level, the high-energy scattering can be described in terms of Wilson lines. In 
addition, from the analysis of quark loops I find that the argument of coupling constant in the BK 
equation is determined by the size of the parent dipole rather than by the size of produced dipoles. 
These results are to be supported by future calculation of gluon loops. 
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I. INTRODUCTION 

At high energies the particles move very fast along straight lines, hence they can be described by Wilson lines U r, {x\_) 
- gauge factors ordered along straight-line classical trajectory of the particle moving with rapidity r\ at the transverse 
impact parameter xj_ (for a review, seeQJ). For deep inelastic scattering, the propagation of a quark-antiquark pair 
moving along straight lines and separated by a distance in the transverse direction can be approximated by the 
color dipole U(x±)W (y±) - two Wilson lines ordered along the direction collinear to quarks' velocity. The structure 
\Q • function of a hadron is then proportional to a matrix element of the color dipole operator 

5: i 

Jl' U r >(x ± ,y ± ) = l-—Tr{U*>(x ± )U^(y ± )} (1) 

switched between the target's states (N c — 3 for QCD). Approximately, the gluon parton density is 

x B G(x Blf i 2 = Q 2 ) ~ (p\lT>(x ± ,0)\p)\*>=Q-' (2) 

• • i 

^ 1 i o 2 

j^j where r\ — In and xb = 2{p-q) IS ^ ne Bjorken variable. 

The small-x behavior of the structure functions is governed by the small-x evolution of color dipoles0,Q. For 

sufficiently small dipoles x\ ~ Q~ 2 so a s (Q) <Cl and we can use pQCD. At high (but not asymptotic) energies we 

can use the leading logarithmic approximation (LLA) where a 8 <C 1, a s hiXB ~ 1. In the LLA, the high-energy 

amplitudes in pQCD are described by the BFKL equation^] leading to the power behavior Fz(xb) ~ x B 12 " ln2 . 
However, the example of DIS from very large nuclei shows that the BFKL equation is not sufficient to describe the 
small-x behavoir of structure functions even in the LLA. Indeed, at sufficiently large atomic number A we get an 
additional parameter a 8 A x ^ ~ 1 which must be taken into account exactly to all orders of the expansion in this 
parameter. The situation is essentially semiclassical: we have a s <C 1 and agF^ ~ 1 where is the strong field of 
the nucleus gluon cloud. Thus we need the LLA in the semiclassical QCD (sQCD): a s <C 1, a s Inxs ~ 1, ot s F^ u ~ 1. 
This situation appears to be general for sufficiently low xb' even for the proton, where we do not have the large 
parameter A to start with, the power behavior of gluon parton density will lead to the huge number of partons in the 
target leading to the state of saturation]^ described by Color Glass Condensate in sQCDj^Q- 
The LLA evolution equation for the color dipoles is non-linear[s|,E|: 

d a N r (r — v) 2 

- U(x, y) = ^ d 2 z 1 J> [U(x, z) + U{y, z) - U{x, y) - U{x, z)U{z, y)\ (3) 
drj 27r z J [x — zj z (z — y) z 
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The first three terms correspond to the linear BFKL evolution and describe the parton emission while the last term is 
responsible for the parton annihilation. For sufficiently high xb the parton emission balances the parton annihilation 
so the partons reach the state of saturation with the characteristic transverse momentum Q s growing with xb as 
e c\nx B _ ^phe ar g um ent of the coupling constant in Eq. © is left undetermined in the LLA, and usually it is set by 
hand to be Q s . Careful analysis of this argument is very important from both theoretical and experimental points 
of view. From the theoretical viewpoint, we need to know whether the coupling constant is determined by the size 
of the original dipole \x — y\ or of the size of the produced dipoles \x — z\ and/or \z — y\ since we may get a very 
different behavior of the solutions of the equation J^J (although first numerical simulations indicate a slow dependence 
of the cross section on the choice of the scale^tj). On the experimental side, the cross section is proportional to some 
power of the coupling constant so the argument determines how big (or how small) is the cross section. The typical 
argument of a s is the characteristic transverse momenta of the process. For high enough energies, they are believed 
to be of order of the saturation scale Q s which is ~ 2 -j- 3 GeV for the LHC collider. Thus, we see that even the 
difference between a(Q s ) and a(2Q s ) can make a huge impact on the cross section. 

The argument of the coupling constant cannot be determined in the LLA so the next-to-leading order (NLO) 
calculation is in order. In the next-to-leading order the non-linear equation @ looks as follows (x,y,z.. are the 
transverse coordinates) 

J- Tr{U x U}} = ^J d2z { a s (x } X ~ 2 ^ 2 _ y)2 + alK NLO (x,y,z))[Tv{U x Ul}TT{U z U f y } - N c Tr{U z Ul}} 

+ a 2 s J d 2 zd 2 z'(K i {x, y, z, z'){U x , Xj\,,XJ z% U}} + K 6 (x, y, z, z'){U x , ul,U z ,,U z , t/j, [#}) (4) 

where -Knlo is the next-to-leading order correction to the dipole kernel and K4 and Kq are the coefficients in front of 
the (tree) four- and six- Wilson line operators with arbitrary white arrangements of color indices. Note that -KTnlo must 
describe the non- forward NLO BFKL contribution found recently in Ref. 11. (The contribution ~ Kq proportional to 
six Wilson-line operators the was obtained in Ref. IT2T) . The calculation of the quark part of the kernel is performed 
in the present paper and the last remaining part of Eq. . the calculation of the gluon part of JCnlo and K± - is 
in progress. 

It should be mentioned that NLO result does not lead automatically to the argument of coupling constant in front 
of the leading term in Eq. 0] In order to get this argument, we can use the renormalon-based approach|l3(: first we 
get the quark part of the running coupling constant coming from the bubble chain of quark loops and then make a 
conjecture that the gluon part of the /3-function will follow that pattern (see the discussion in Refs H . 

As we demonstrate below, the result is that the value of coupling constant is determined by the size of the original 
dipole rather than the size of the produced dipoles: 

±U(x,y) = a "^ x ~y) 2 ) N - J # z } X zn y f_ y)2 Mx,z) +U(y,z) -U{x,y) -U(x,z)U(z,y)] + ... (5) 

The paper is organized as follows. In Sect. 2 I recall the derivation of the BK equation in the leading order in a s . In 
Sect. 3, which is central to the paper, I calculate the quark contribution to the small- a: evolution kernel of Wilson-line 
operators. In Sect. 4 I present the arguments that the coupling constant in the BK equation is determined by the size 
(x — y)j_ of the parent dipole. The light-cone expansion of the quark-loop propagator is performed in the Appendix. 



II. DERIVATION OF THE BK EQUATION 

Before discussion of the small-x evolution of color dipole in the next-to-leading approximation it is instructive to 
recall the derivation of the leading-order (BK) evolution equation. As discussed in the Introduction, the dependence 
of the structure functions on xb comes from the dependence of Wilson-line operators 

lT>(x±) = PexpUg f dupl A^iupV +x±)\, p" = pi + e~ ri p 2 (6) 
J —00 

on the slope of the supporting line. Here p\ and p 2 are the light-like vectors such that q = pi — xbP2 and p = 
2 

P2 + — pi where p is the momentum of the target and m is the mass. Throughout the paper, we use the Sudakov 
variables p — ap\ + (3p 2 + p± and the notations x m = x^p^ and x* = x^p 2 related to the light-cone coordinates: 

To find the evolution of the color dipole QJ with respect to the slope of the Wilson lines in the leading log 
approximation, we consider the matrix element of the color dipole between (arbitrary) target states and integrate 
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FIG. 1: Leading-order diagrams for the small- a; evolution of color dipole. 



(d) 



over the gluons with rapidities rji > 77 > ?/ 2 = Vi ~ ^rj leaving the gluons with rj < 772 as the background field (to be 
integrated over later). In the frame of gluons with 77 ~ 771 the fields with 7/ < 7/2 shrink to a pancake and we obtain 
the four diagrams shown in Fig. ^ 

The (background-Feynman) gluon propagator in a shock- wave external field has the form[iJ^| 



(A^(x)Al(y)) = 9(x*y*)(x\ g ff \y) 
^ l(P 2 + it) 



{x*)0(—y*) / da 



Aa 2 



{ Xl _\e- l — x * 2a 9llv U+ -{id^U - p 2fl id u U) - ^ djU 



-x*)0(y*) / da 



> e ia(x-y). 



Aa 2 



(x±\e 



2ag^ - -{id^p 2v -p 2ll id v tf) - ^E^Q 2 ^ 
s as 2 



(7) 

^ v *\y±) ab 



where d\ = —did 1 . Hereafter use Schwinger's notations (x\F(p)\y) = fctp F(p) e ip '^ x v > and (x±\F(p±)\y) = 
J dp e lt -P' x ^ v ' ,± ( the scalar product of the four-dimensional vectors in our notations is x-y = | (x*y, +X*y») — (x, y)±) 
We obtain 



du I dv A a ,{up" n + x x )A b .{vp r > 1 + y±)} 



Fig[TJi 



, da 

= - 2a s / — [x± 



2 TT ab 



diU 1 



-12 ol p 2 l +a 2 e 2? ' 1 s p 1 ^ + a 2 e~ 2r >i s 



\v±) 



(8) 



Formally, the integral over a diverges at the lower limit, but since we integrate over the rapidities 77 > 772 we get in 
the LLA 



00 rO 



du / dv A a m (up v + x±)A b , (vp v + y±_)\ m 
J -co JFigia 

= -2a s Ar,(x x \±-d 2 ± U ab ±-\y ± ) = - 2a s A V f d 2 z ± (x ± \^\z ± )(2U z - U x - U y ) ab (z x \^-\y±) 
Pi Pi J P \ P\ 



(9) 



and therefore 



2tt 2 



Ary {t a U x ®t b Ul) m d 2 z 



(x - z,y - z)_ 
(x - z)\(y - z) 



2 vu?-u?-u?)<* 



(10) 



The contribution of the diagram in Fig. \\]p is obtained from Eq. I|10|) by the replacement t a U x ® i b V\ — > U x t b (g> Ult a , 
x «-> y and the two remaining diagrams are obtained from Eq. El by taking y = x (Fig. ^p) and x = y (Fig. CQ})- 
Finally, one obtains 



{ u *® u y}v lg m = -^{ ta u*®t b ut + u x t b ®uln^Jd 2 z {x z > y z) 



■ '^{'"i , r; r I (up - uT) ab + ^P-{u x ® t b ult a y* J £J± 2 (c/f - ^)° b (11) 



d 2 z 



2 -U? -Up)* 



For the color dipole Q one easily gets the BK equation 
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III. QUARK CONTRIBUTION TO THE NLO BK KERNEL 
A. Quark loop in the momentum representation 

There are two types of quark contribution in the NLO: with quarks in the loop interacting with the shock wave 
(see Fig. or without ( Fig. |2p). (In principle, there could have been the contribution coming from the quark loop 
which lies entirely in the shock wave, but we will demonstrate below that it vanishes). 

The quark propagator in a shock- wave background has the form [j| : 



tp{x)tp(y) = 0(x*y*)(x| 



p- + ie 



\y) 



da 
2a 2 s 
da 



2a 2 1 



e -i<*(*-y). (a^Kopi + p ± )e-^ x *p 2 Ue^y' {api +P±)\v±) 
-ia(x-„). ( x ±\(afa. +p±)e- l ^ x *p 2 U^e^ y *(ap 1 + p±)\y±) 



(12) 



ipt a M(x)ip(y)^ 1 i/;(y) 
da f 1 dv 



Multiplying two propagators one gets at X* > 0, < 

^ i^-y). trTr^l^xe-^^e^^^^^l^e^.^teilL^^I^) (13) 

Iq lbira 3 J v z v z 

where tr stands for the trace over spinor indices. Therefore the quark-loop contribution to the gluon propagator is 
(At{x)A b M) 



— 2a a n 



'da 



Tr{t a U z t b ul,} I dv 



d 2 zd 2 z'J d 2 kid 2 k[d 2 k 2 d 2 k' 2 e i(ki,^)±+i(k[^-z')±-i(k2,y- Z )±-i(k' 2 ,y-z') ± 

(A* , k 2 ) (hi , k' 2 ) + {h , k[) (k 2 , k' 2 ) - (h , k' 2 ){k[ , hi) 



[(fei + k'^vv - kjv - k'^v][(k2 + k' 2 ) 2 vv - k\v - k' 2 v 



. (fci+fci)j , 



(14) 



and one obtains the contribution of the diagram in Fig. [2K in the form 



( / du f dv A a ,{u Pl + x±)A b m (vpi +y±)\ 
^Jo J -00 )Fi t 



ia 2 s n f 



'da 



— jdzdz' / d 2 k 1 d 2 k[d 2 k 2 d 2 k' 2 e ^.x-z)^+ l (k' 1 ,x-z')^- l (k 2 .y-z)^- l (k' 2 . V -z') 1 _ 

(fci,fc 2 )(M,^) + (k u k[)(k 2 , k' 2 ) - (fr, to) 



Tr{t« Uz t b ui,} Cdv <Mm 

Jo (ki + k[ 



) 2 {k 2 + k' 2 ) 2 {k 2 v + k' 2 v){k 2 v + k' 2 v) 



(15) 



where n/ is a number of light quarks (n/ = 3 for the momenta Q s ~ 1 -j- 2 GeV) and Tr stands for the trace over 
color indices. The variable v is the fraction of the gluon's momentum a carried by the quark. 

To calculate this diagram we use the dimensional regularization and change the dimension of the transverse space 
to d = 2 — 2e. The calculation yields 



4a? 



-71/ 



A V Ti{t a U x t b Uh d d pd d qd d q' 



e i(p,x)x-i(p-q-q' ,y)i 

p 2 {p — q — q') 2 



dvdu 



-(q + Q') 



'\2 



uuT(e)jjL 



2 c 



(P 2 vv + Q 2 uuY 



r(l + e) M 2 



[P 2 vv + Q 2 uu) 



l+e 



{p 2 [vv(q,q') - uuQ 2 + 2uuvv(q 2 + q' 2 )} - 2uuvv{P, q){P, q') 



+ uu(l - 2u)[vv(q, q')(P, q + q') + vq 2 (P, q') + vq ,2 {P, q)} + u 2 u 2 Q 2 (q + q') 2 } 



Tr{t a U{q)t b U\q')} 



(16) 



(a) (b) (c) 

FIG. 2: Quark-loop contribution to the gluon propagator in a shock-wave background . 
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FIG. 3: 



where P = p — (q + q')u, Q 2 = q\v + q'^v and we use the notation J ct d p± = J d ^ d . 
The contribution of two diagrams in Fig. |2t>,c is 

2n f a 2 s A V tr{t a U x t b Ul}^B(2 - e,2 - e )(x x \±{^, dlU ab }±\y x ) 

I" P± L (Pi) J P± 



(17) 



where B(a, b) — r(a)r(fo)/r(a + b). The contribution of diagrams in Fig. 131 is obtained from the sum of Eq. (|16fl an d 
(11711 by the x <-> y replacement in the coefficient in front of tr{t a U x t b U^} and the contribution of the diagram in Fig. 
0] by taking y — > x in this coefficient and changing the sign. Similarly, the diagram in Fig. is obtained by taking 
x — > y. The sum of all diagrams has the form 



^n f Ar,Tv{t a U x t b Un 



f <t d p S d q efV , i(PtX)± _ e i(p, v ),w e -i( P -,- g ' lJ; ) 1 _ e ~i{v- q - q ',y)±) dudv 

p ip-q-Q) Jo 
r(i + e) M 2 - d 



,. 2 4^r(e)/z 2e 



(P 2 vv 



i ( „ 

— \P 2 [vv{q, q') - uuQ 2 + 2uuvv{q 2 + q' )] - 2uuvv(P, q)(P, q') 



+ uu(l — 



(1 - 2u)[w(«, 9 + <?') + ^(P, 9) + W, g)] + u 2 u 2 Q 2 {q + g') 2 }] Tr{t a U(q ± )t b UHq' ± )} 



2^B(2 e, 2 e) [ 2 ( Xx \-L{^l, fl^}-^) - 

(y^{^dlU« b }-L\y ± )} +±-[ 2 (x ± \^dlU« b -L\y ± ) - {x± \-L d lU* b ^\ X± ) ( y± \-LdlU« b -L\y_ 



-\y ± ) - (x ± \^dlU ab ^-\x ± ) - {y^dlW 



(18) 






FIG. 4: 



FIG. 5: 




FIG. 6: Quark loop inside the shock wave . 



where the last term ~ J- is a counterterm calculated in the Appendix. 



B. Quark loop inside the shock wave 



Let us consider the diagram in Fig. 0with the quark loop inside the shock wave of width A = Az» ~ -^e 7,2 where 
m 2 is some characteristic mass scale of order of Q 2 . From the form of the perturbative propagator 



p z + it 8ttx* Jq a 

we see that the characteristic transverse scale inside the shock wave is 



A e~ m 

{z-z')\ — m~ <rn~ 2 (19) 



as a 



and therefore the contribution of the diagram in Fig. reduces to the contribution of some operator local in the 
transverse space. By dimensional arguments, this local operator must have the same twist as the operator describing 
the interaction of the gluon with the shock wave at the tree level. In the leading order in a s , the vertex of interaction 
of gluon with the shock-wave field is proportional to 

d\U = - ig[DG] + 2g 2 [GG] (20) 

where 

[DG] = J du [oopi,upi] x D l Gi.(upi + xj_)[upi,-oopi] x (21) 

[GG] = J dudv 8(u~v) [oop 1 , upi] x G mi (upi + xj_)[upi, vpi] x G, l (vpi + x±)[vpi, -oopi] x 

These operators have twist 2 so a possible local operator describing the gluon interaction with the shock wave at the 
one-loop level must also be of twist 2. To find this local operator, we consider (the quark loop contribution to) the 
color dipole tr{U x U^} at small at (x — y)\ — » and compare the expansion of the contribution of the diagrams in 
Fig. [21 to the exact calculation of the light-cone expansion of tv{U x Ul} in QCD (up to twist-2 level), see Fig. 

The first step is the light-cone expansion of the sum of the diagrams in the in Figs, yj- |S| in the shock-wave 
background. The light-cone expansion of Eq. (|18[) at x± — * y± starts with terms quadratic in q(q'). They lead to the 






.... 





--0+ 



FIG. 7: A possible local contibution coming from the quark loop inside the shock wave 



operators d 2 U x and diU x diU\ (in the leading order we do not distinguish between U x and U y ): 



Ti{diU x diUl} x ^ ^n f Ar)U : 



jab 



fi 2 - d B(l- e ,l-e)(x ± \{[3 



i 1Jt PiP 3 \r(i + £) 



\y A _)Tr{1?d i U x t b d j U} 1 } 



lp_lJ Je -Pi 

a: . 1 ,,„„ . slYl 



v 2 

7T 
,,2e 



-3 



1 



r(-2e) 1 r(-e) 



e S(2 £ ' 2 e) (Ai)"2e i2e(A5_)-« 



2(1 + e) 



, 2e AiAj 

tin H — 



(l + £ ) Ai J(Ai) 



r(-2 £ ) 



TridiUxdjUl} 



(22) 



The expression (1221) should be compared to the light-cone expansion of the quark-loop part of the gluon propagator 
in an external field (see Fig. 01 performed in the Appendix. A typical term of the light-cone expansion has the form: 



AZ(z,,xx)A b m (y.,yx) 



ig 3 f°°da/4iy 



(23) 



[ x * 2 1 

ig I d-z* -\(x - z)% + (z - y)%] z*] x D l Gi,(z*, xx)[z*, z'*\ x {z*, y*]x) ab 



+ g62^ '<£z*J 'd-si ({x*,z4 x G:(z*,x ± ){z*,z:} x G. l (z:,x ± M,y4 x ) ab {^{(x - z)l + (z' - y)%] + (z - z')%] 

In our "external" field the characteristic distances are of the order of width of the shock wave: z*,z* ~ 

e m y's/m 2 . As we shall see below, the characteristic distances x* and y* are ~ e m yjs/m 2 so we can neglect z* and 
z' t in comparison to x* and/or y*. The formula (|69[) simplifies to 



1d7t z ./n a J \as 



)( 



47TA, 



1-6 A i 



(24) 



1 Z^* 2 

«5- / d-z* [xl + (-y)l]([x*,z*]xD*Gi.(z*,xx)[z*,zi] x [z*,y*] x ) c 



+ g 2 [ d~z* [ d-zi {[x*,z4 x G.\z*,xx)[z*,zl} x G mi {zi,xx)K,y*} x ) ab {-[xl + (-y)l] + (z - z')l) 
Jy. s Jy, s Le J 

A very important observation is that the contributions proportional to 

g 4 n f [dz* f dz 1 , 6{z - z') (z - z') e G.i{z*)G.i{zi) (25) 
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present in the individual diagrams in Fig. ^] cancel their sum. If it were not true, there would be an addtional 
contribution to the gluon propagator JJJ at the g 4 level coming from the small-size (large- momenta) quark loop. 
Indeed, the calculations of Feynman diagrams with the propagators (0 and l|12[) imply that we first take limit 
z *: z '* — * and limit d± — > 2 afterwards. With such order of limits, the contribution (|74|l vanishes. However, the 
proper order of these limits is to first take d± — > 2 (which will give finite expressions after adding the counterterms) 
and then try to impose condition that the external field is very narrow by taking the limit z*, z' r — > 0. In this case, 
Eq. I|74[l reduces to g 4 n f [GG] . The non-commutativity of these limits would mean that the contribution \ [GG] \ 
should be added to the gluon propagator (JJJ to restore the correct result. Fortunately, the terms ~ <|74l) cancel which 
means that there are no additional contributions to the gluon propagator coming from the quark loop inside the shock 
wave (= quark loop with large momenta). 

Since there is no external field outside the shock wave, after cancellation of the terms ~ (z — z') e we see that at 
> Eq. (|24|) vanishes, and at x* > > one can extend the limits of integration in the gluon operators to 
±oo and obtain 



W 

167T' 



■B{2 



of *?(-) 

.In a 6 \as/ 



las 
4ttA, 



(26) 



■ig-K + (-y)l][DG\ 



ab 



+ g ^GGr-[xl + (~y)t 



167T 



:B(2-e,2-e) 



' d~a f Ai 
a 3 \as 



las 
47rA» 



[< + (-y)l] diu a x 



ab 



The light-cone expansion of gluon propagator contains only Wilson lines and their derivatives as should be expected 
after cancellation of the "contaminating" terms iJUJ- 

We have demonstrated in the Appendix that the light-cone expansion of the quark-loop contribution to the gluon 
propagator coincides with Eq. lf2~2")l as should be expected once we established the commutativity of the limits d± — > 2 
and z*(z'^) — > 0. 



C. Quark loop in the coordinate representation 



To calculate the integrals over momenta in Eq. I|18|) it is convenient to subtract (and add) Tr{t a U z t b U].} from 
Tr{t a U z t b ul,}: 

Tr{t a U z t b ul} = Tr{t a U z t b ul - t a U z t b Ul) + Tr{t a U z t b Ul} (27) 

Let us start with the last term in the r.h.s. of Eq. 127|) . In the momentum representation, this term corresponds to 
Tr{t a U(q)t b C/t (<?')} 4tt 2 (%') / dze l ^' z ^Tr{t a U z t b W z } = 2n 2 S(q')U ab (q) so we get 



. , p i(p,x)±-i(p-q,y)± 

-2a 2 s n f ^T 1 t a U x ®t b Ul^- U d P <t d - 



p2(p _ q y U»>( q )ldvdu 



uuq 2 T(e) 



[(p — qu) 2 vv + q 2 vuu] 



q 2 uuvT(l + e) 



[ip — qu) 2 vv + q 2 vuu 

■2 



— {(p - qu) 2 (-l + 2v) + uuq 2 } 



2^n f ^t a U x ®t b Ul d- d pct d q 37—^ U ab (q) q 2 , _ 2 ) ^ e B(2 - e, 2 - e) 



p 2 (p — q) 2 y ' 1J * (q 2 1 'n 2 y y ' ' ^ ^ 

where we have used integration by parts to transform the second term in the l.h.s of this equation. Alternatively, this 
result can be obtained directly from Eq. Ijl5(l after the substitution H27JI . 
For future use we need to rewrite it in Schwinger's representation: 

■2 



2^f- d n f Ar 1 {t a U x ®t h Ul}{x\\ ( F(e) « 2 "«^ 1 



p 2 v(-a 2 ) e 



diU ab )-\ y )B{2-e,2-e) 



V 



(29) 



The contribution coming from the first term in Eq. (|27|l is UV-finite. To calculate it in coordinate representation it 
is convenient to return back to the original expression 



5a„ 



n f Ar]t a U x ®t b Ul j d- 2 k 1 <t 2 k' 1 d 2 k 2 d 2 k' 2 e ^^-z)x+i{K,w-z')x-KKy-z)x-i{k' 2 ,y-z') ± 

Tr{t a U z t b U\ - t a U z t b U^} t du fc2)(fci ' k>2) + ( k ^ k 'i)(k2,k' 2 ) - (h, k' 2 )(k[, k 2 ) 



o (fci + k' 1 ) 2 (k 2 + k' 2 ) 2 (kju + k[ u){klu + k' 2 u) 
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and use the formulas 



****** (^)^-)+^ 2 ) 



(fci + k 2 ) 2 (kfu + k\u) Att 2 (xi - x 2 ) 2 



**h**k 2 



kuh 



li«2j - * <-> J 



i(k 1 ,x 1 )+i(k 2 ,x 2 ) 



1- I + 0(d-2) 

xfu + x 2 ui 

x u x 2j -i<-> j 

i7r 2 (xi — x 2 ) 2 (x\u + x\u) 



+ 0(d-2) 



(fci + k 2 ) 2 (k 2 u + k 2 u) 
After some algebra, one obtains: 

(x,x') (x,x')(r ! r') + (x,r)(x , y , )-(^^')(^'^) 



(30) 



U x ®Ul = - ^jA V {t a U x (g> iW} / d 2 zdV / Tr{< a [/ 2 i b [/j, - i a C/ 2 i b [/]} 



1 



1 - 



X 2 u + X' 2 u Y 2 u + Y' 2 i 



(X 2 u + X' 2 u){Y 2 u + Y' 2 i 



(31) 



where X = x — z, X' = x — z' ,Y = y — z, y' = y — z'. Note that the singularity at z' = z is integrable. 

Performing the integation over u and adding the {/^-divergent term (|29|l one obtains the total contribution of the 
diagram in Fig. [2i in the form: 



2tt 4 



Ar){t a U x ®t b Uh / d 2 zd 2 z' Ti{t a U z t b ul - t a U z t b U\} 



(z - z'Y 



x U 



(x,x') x 2 (y, y') ^y 2 t (i,i')(y,y') + (i,r)(ir) - (i,y')(^7) ^ y 



111 



J^"2 J£> 2 j£i 2 _ y 2 y 2 



2^n 2e B(2 - e, 2 - e)A?7{t £/ x ® t^ttf }(icx 



1 / F(e) 



piV(-5i) E 



3 2 (7 



X /2y2 _ Y' 2 X 2 
ab I 

—\y±) 
pi 



in 



y/2 X 2 



Sum of the UV-divergent contributions takes the form 

< 2e 1 r r(e 



The counterterm is calculated in the Appendix (we use the MS scheme): 

1 , , 1 „,„„ % _1 

p\ 

Adding the counterterm, one gets after some algebra 

n f ^-A V t a U x ® t h Ul{x x \±- [{ In ^, a 2 [/} - f In -i^U 2 7 ^ -i- ' 

'g? 2 z r (x — y) 2 



-2n f alAr,t a U x ®t b Ul — { Xl \^dlU ab ^\ yi _) 



(32) 



(33) 



(34) 



P± 



1 

X 2 



»t£Wu,**uij 8n2 L x2y2 
i 



{- 



[ln(.x - y) V + 



[lnF V + -] - — 2 [InAV + -]}(2L/ Z - C/ x - C/,) Qb 



(35) 



The calculation is simplified if one notes that &]U Z in the l.h.s. can be repaced by d±(U z — \U X — \U y ). Our final 
result for the sum of the diagrams in Fig. [2] has the form 

(U x <g> Ul)^ = ^A V n f {t a U x ® t b Ul) jdz^-J dz' Tr col { - t a U z t b Ut) j-L^ 

j jx,x') 1 a 2 (y,y f ) , y 2 (a, a')(f, rQ + (a, y)(x'y') - (a, f / )(a / , y) , A ,2 y 2 ] 

x si In — ft tt In — o ~~ \~ o o In ■ 



A /2 r 2 - r /2 A 2 y /2 A 2 
Y 2 {^x^ Mx " y)V + 1 ~ A^ [lny V + 1 " ^t lnX ^ 2 + ^}( 2 ^ " ^ " 



y2 _ JT' 2 Y' 2 ^2 _ y/2 y/2 
^2 



(36) 
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As we mentioned above, the contribution of diagrams in Fig. |3]is obtained from Eq. (|31)|) by replacement t a U x ®t b U^ y — > 
U x t b <S U^t a and the contribution of the diagram in Fig. 0] can be obtained from Eq. (|36[1 by taking y — x in the 
integrand (and changing the sign) 

(U x ® ^) F . gE | = ^fAry n f {l»Uj> ® ^d? ± | dz' Tr{t a U,t h Ul - t a U z t b U\) 



x 2 -x' 



Similarly, the contribution of the diagram in Fig. [SJis obtained by the replacement 



a; -> y: 



{u x ® 



x 1- 



FigEl 

(y,y) 



A?7 n/^®* 6 */***} / da! 



- / d»' Tr col {t a U z t b U f z , - t a U z t b Ul) 



(z - z'Y 



Y 2 



Summing the contributions of the diagrams in Fig. [2|-[S|and taking Tr over the color indices, one obtains 

2 f 

Tt{U x Ul) = ^Ar)n f Tr{t a U x t b U^}J d 2 z 



- [ d 2 z' Tr{t a U z t b ul - t a U z t b Ut) . 1 A : 
ttJ z (z- z'Y y 



X' 2 Y 2 + Y' 2 X 2 -(x-y) 2 (z-z')\ X' 2 Y 2 



2(X' 2 Y 2 - Y' 2 X 2 ) 



In- 



Y' 2 X 2 



(x-y) 2 „ , ,22 5 n X 2 -Y 2 , X 2 

X 2y2 W^-») A 4 +3]+ X 2y2 ln y2 



Let us present the total result for the sum of the leading order BK equation and the quark NLO correction 
£jT±{U s U}} = ^Jd 2 z [Tr{U x Ut}Tr{U z U%} - N c Tr{U s l4}] 



(x - yf 



X 2 Y 2 







£5!»£m*-v)V + - 

6tt V yi ' 3 



Utt x 2 y 2 n r 2 



4n/'n{* Cxt 6 Dj} / d 2 zd 2 z' Tr{t a U z t b ul, - t a U z t b Ul) 7^-^ 



{ 



x 1- 



X' 2 Y 2 + Y' 2 X 2 



(x-yY(z-z') 



l\2 



2(X' 2 Y 2 - Y' 2 X 2 ) 



In 



X' 2 Y 2 



Y' 2 X 2 f 



(37) 



(38) 



(39) 



(40) 



We see the first term proportional to ln(...) 2 /i 2 (we will call it the "UV" term) has the same structure as the 
zero-order contribution In the next Section we will use it to determine the argument of the running coupling 
constant in Eq. ©. 



D. Comparison to NLO BFKL 

To compare with the NLO BFKL equation we need to linearize Eq. (I40|l which gives 

d 2 z iU{x,z) + U(z,y) — U(x, y)] 



U{x,y) = ^ 
dr] 27T Z 



(x - yf 



X 2 Y 2 



()7T 



a s n f X 2 - Y 2 
~6tt X 2 Y 2 



In- 



X 2 
Y 2 



a l n f 
4N r ir 4 



2 , U(z,z') r 

zdz yT^yX 1 



X' 2 Y 2 + Y' 2 X 2 -{x- y) 2 (z - z'f X' 2 Y 2 ' 



2(X' 2 Y 2 - Y' 2 X 2 ) 



In 



Y' 2 X 2 - 



(41) 



This should be compared to the quark part of the non-forward BFKL kernel |17| but the Fourier transformation from 
the momentum space to the dipole-type representation appears to be rather difficult. 
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To simplify the comparison, let us consider the case of forward scattering and write down the Mellin representation 
oiU(x,y) 



Wlv-u) = \dv (x-y) 2 ~<W u , 1=\+iu 



(42) 



where we have displayed the dependence on the rapidity 77 explicitly. Using the integrals (x(t) = [ — V'(t) — ^(1 — 7) + 
2V(1)]) 



2 _ (x-y) 2 r 



d 2 z 



(Pi 



X 2 Y 2 

\nX 2 /Y 2 \nX 2 /Y 2 



(43) 



Y 2 

V 



(z - z') 4 - 27 
we obtain 



1 - 



X 2 
X' 2 Y 2 



{X 2, + y2 1 y) 2 T) = n{x _ yfn^'il _ 7 ) _ ^'( 7 ) _ x 2( 7) + _ x(7 ) } 

7 



Y' 2 X 2 



(x - y) 2 (z - z' f 



2{X' 2 Y 2 -Y' 2 X 2 ) 



In 



X' 2 Y 2 
Y' 2 X 2 



7T COS7T7 2 + 37(1 —7') 

sin 2 7r 7 (l-4 7 2 )(3-2 7 ) 



(x - y) 



2~ 



5* 



OsNcr a s 2 2 

7T 07T 



(44) 



(x(7) ~ T^"/{[^'(7) - 1^(1 - 7) + X 2 (7) - -X(7)] + yX(7) + 



37r 2 cos7T7 2 + 37(1 — 7) 
N 2 sin 2 7r 7 (l-47 2 )(3-27) 



where 7 = i + iv. This expression should be be compared to the NLO BFKL result 18]. Unfortunately, there is no 
explicit expression for the coordinate-space NLO BFKL kernel yet. However, the last two terms in braces in r.h.s. 
of this Equation coincide with the expression for the n/ part of the eigenvalue 5(7) of Ref. 18; The first term in 
braces should correspond to the quark part of /3-function contribution to the eigenvalue £(7). We expect to study the 
relation to NLO BFKL in detail after completing the calculation of the gluon loop. 



IV. BUBBLE CHAIN AND THE ARGUMENT OF COUPLING CONSTANT. 

To get an argument of coupling constant we can trace the quark part of the /3-function (proportional to n/). In 
the leading log approximation the quark part of the /3-function comes from the bubble chain of quark loops in the 
shock- wave background (cf. Ref Il9|) . We can either have no intersection of quark loop with the shock wave (see Fig. 
[SJ or we may have one of the loops in the shock-wave background, (see Fig. |§J) 

It is easy to see that the sum of these diagrams yields 

^Tr{U x Ul} = 2a s Tr{t a U x t b Ul} Jd 2 pd 2 q [ e i( - p ' x ^ - e *<P>v)±][ e -*<P-**)± - e -<(p-9.w)x] 

x ^fl-^ln^iCT^) 1 (45) 

P 2 (l + fjln^)V (P-ff) a (l + £k<i£j*) 

where we have left only the UV part l|29|) of the quark loop. (In principle, one should also include the dressing of the 
UV-finite 1/N C term in Eq. (|41|) by bubble chain, but I think that it is a separate contribution which has nothing 
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to do with the argument of the BK equation). Replacing the quark part of the f3- function — ^n/ln^ by the total 
contribution fjMn ^ (where b — jjq we get 



Tr{U x Uh = 2Tr{t a U x t»U}} 



x / d 2 pd 2 q [e^ x ^ - e ^y^}{e~ l(p -«< x ^ - e - i ^-^]^^a- 1 (q 2 )dlU a \q) ^- q P (46) 

J ' v (p-qr 



To go to the coordinate space, let us expand the coupling constants in Eq. (|4(j[l in powers of a s = a s (/i 2 ), i.e. 
return back to Eq. (|43|l with §Jn/ — > — &§^. In the first order we get the UV part of the NLO BK equation (flOf 



(47) 



2tt 2 



d 2 z [Tr{U x Ut}Tr{U z Uh - N e Tr{U x lft} 



(x - y) 2 



a s X 2 - Y 2 X 
In ■ 



2 -, 



In this section, we perform the calculations in the leading log approximation 

-,2 



, In 



El 

/< 2 



l,a s < 1 



(48) 



(hence we omit the constant term (~ | ) from the Eq. iJlDJ) 
In the second order in the expansion we obtain 



2a s / d 2 pd 2 q e <(i'.*)j.-<(p-«.»)x 



47r g 2 



4F^^ - (f-9) 2 (l-^ln^) 

fe 2 2 In 2 4 1 In 4 In 4 1 In 2 4 hi 4 / u 2 \ 1 1 / u 2 \ In 4 



8tt 2 



,2 "^i^ 



-a 2 f/- 



P A p* P A 



In- 



-9 2 



d 2 £/ )-=■ - — « f In 



p2 ^2 



-9 2 



Again, it is convenient to replace d^U z by d\U z where U z = U z — — 
Using the formulas 



(x\ - In 2 



p2 p2 p2 



[/In 2 



1 1 



^ — 9 2 P 2 ID 2 — f' 1 "' 



u 2 - In ^ 

pZ pZ 



In ^ - u 2 

pZ pZ 



—dj_ ) p 2 p 2 



2 2 



2 



»> + M-i*-M[( ta :V). + ( h V). 



1 



l47T 2 X 2 



X 2 



In X V In Y V - In — In A V - In 2 A V 

y 2 



47r 2 r 2 



lnXVlnr 2 ^ 2 +ln^lnA 2 /i 2 - In 2 A 2 /x 2 ] }((7 2 - hj x - Uj y ) 
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and 



(x\2 



Pi In 



2 ^ 



^U^r + 2 



Pi In K _pjln 



V —U- 



P 



P 



Pi - Pi In 2 

9 9 



9 Pi 
Z 9 
P~ 

A 2 

x 2 r 2 



In- 



-9 : 



Pi In 



ln 2 AV 



we obtain 
d 

drj 



\v) 



2Aj 
ttA 2 



Tr^C/t}^) = ^(^y l^ z[ TT{U x UnTT{U z Ul}-N c TT{U x Ul}} 



dz 



1 



;t4 



i 



dr 



An 2 y ~- 2 2 
ri Y 2 ~(Y,r)Y t 



In- 

r 1 1 

X 2 + Y 2 



-d- 



-U 



ln 2 AV 



r 2 (Y ~rf 



(49) 



(50) 



ln 2 (x - y) V + ± In fj(ln* V + ln(z - y) V) - ^ In fJ(lnF V + ln(x y) V) 

We have omitted the contribution of the last integral in r.h.s. of Eq. H49|l since it is negligible in the limits X»7, 
Y 3> X and X,Y ^> x — y, and therefore can be dropped in the leading log approximation Q48[l . Adding the first-order 
contribution (first line in the Eq. ijiDjl. we get 



Tr{C/^t } = ^|d 2 ^ r&{^}lt{^}-JV c T>{^}]{ 



drj 

+ (^)V (X _ 9)V 



ba a 1 



X 2 



In 

4tt x 2 y 2 



6a, 



1 + ^ ln(a: - y) V + lnX V 

47T V ^ 47T 



(x ~ yf \ ba s 22 



ba s 1 X 2 - 
— - — In 

4tt Y 2 Y 2 



;i + £ M *- y) V + £lnYV]} 
Our guess for the argument of the coupling constant in all orders in \np 2 //i 2 is [23 



6a „ 



-Tr{£/ x [/t} 



"a ((a: - y) 2 ) 



-(z-y) 2 , 1 /a,(X 2 ) 



d 2 z [Tt{U x UI}Tt{U z UI} - N D Tr{UJJ$}] 



(51) 



(52) 



X 2 Y 2 X 2 \a s (Y 2 ) 



-1 



1 /« S (F 2 ) 



F 2 Va s (X 2 



- 1 



We see now that the argument of the coupling constant in the BK equation is size of the original dipole (x — y) 2 as 
it was advertised in Eq. JSJ. 



V. CONCLUSIONS AND OUTLOOK 



First, there are no new operators at the one-loop level - just as at the tree level, the high-energy scattering can 
be described in terms of Wilson lines. The fact that there are no new operators at the one - loop level is rather 
remarkable. In the case of the usual light-cone operator expansion this is not true - for example, if we have the 
operator 

[GG] = Jdudv0(u-v) [oop 1 ,up 1 ] x G, i (upi + x±)[upi ) vp 1 ] x G*(vp l +x±)[vp 1 ,-oop 1 ] x (53) 

in the leading order, one should expect the operator 

a s Jdudv 9(u- v)ln(u-v) [oopi, upi] x G.i(upi + x±)[up 1 , vpi] x G*(vpi + x±)[vpi, —oopx] x (54) 

in the NLO (in general, any new loop brings an additional factor a s ln(w — v)). This does not happen here, and in 
addition the operator [GG] appears only in the combination —i[DG] + [GG] = d±U, exactly as at the tree level. I 
have checked this by the explicit calculation of the quark-loop contribution and expect to confirm it by the calculation 
of the gluon loop. 
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FIG. 10: "Target" contribution to the gluon propagator in the external field . 

Second conclusion of the paper is that the argument of the coupling constant in the BK equation (obtained from 
the renormalon-based arguments) appears to be the size of the parent dipole rather than the size of produced dipoles. 
I have obtained the result for the argument of the coupling constant in the non-linear evolution of dipoles using the 
quark part of the /^-function. It is necessary to confirm this result by calculating the diagrams with gluon loops. 
Also, it would be extremely interesting to check how (and if) this argument of the coupling constant arises from the 
correlation function of the original dipole and the "diamond" high-energy effective action 22] formulated in terms of 
the (renorm-invariant) Wilson lines. The study is in progress. 
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VI. APPENDIX: LIGHT-CONE EXPANSION OF THE QUARK-LOOP CONTRIBUTION TO GLUON 

PROPAGATOR IN THE BACKGROUND FIELD 

The expression <|22[) should be compared to the light-cone expansion of the quark-loop part of the gluon propagator 
in an external field. The quark-loop contribution to the propagator of a gluon in the external field has the form: 

WW =J dx i dy > n p2g ^ - JL- - a- i^)-Tr { ^(^|i|,o^(,^i^)}(,i PV> + 2 ] Gvm + Qvm \ v r 

(55) 

where 

P . = 0.„ = -D ! G V = ^D'G,. (56) 
a as 

An additional term in the gluon propagator is due to the fact that the external gluon field of the target satisfies the 
Yang-Mills equation with a source D^G^ = — gip^^t a ip. From the viewpoint of Feynman diagrams in the bF gauge, 
this term comes from the diagrams with the quark insertions shown in Fig. I1UI (in the light-like gauge this term arises 
automatically, see Ref. |9|) . For the contribution ~ the quark propagator reduces to 

7p {a p + a k ){fi p + p k )s-(p + k)l Pl Pl (a p -a k )(l3 p -(3 k )s-(p-k)l ^ [0<) 
As explained in Ref. Q at at <C a p one can shift the contour of integration over (i p away from the pole in the 

2 

denominators in the above equation. After that f3 p ~ so one can neglect the terms proportional to transverse 
momenta in the denominator and in the numerator. One obtains 

— 7m A At a t h - t b t a fa Mm - -*f abc t c ^ (58) 
as a 

which corresponds to the vertex of the insertion of O^, operator. 

We need to expand the Eq. I|55|) near the near the light cone x — > y and compare it to the light-cone expansion of 
the same propagator in the shock- wave background (refvesvkladlikone). The technique for the light-cone expansion 
of propagators in external fields was developed in Refs. I20ll2l|) . The expansion of the tree- level quark propagator has 
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the formpH: 



1 . AT(2 - e) r(i - e) 



— J du |it AVG(x u ) + ucrG(x u ) A" — liuu A'D x G\a 



( X '\ — \y') = 

{ l f> ]yj 2tt 2 (-A 2 ) 2 - £ ' 16tt 2 (-A 2 ) 
+ 4A^ dv uvG£(x u )G A i:(x v )} 

+ Y^p^p^ du - uu)D x G Xp Y + l -uu{l - 2u)DD x G XA - uu(l - 2u)Y[G£(x u ), G p5 (x„)] 
^uue^uxD^D^^ + i^{[G^{x u ),G}{x u )] + [G M (x u ),Gf(x v )])w 6 



2 — ^-/-'^^ s, i /u - - 2 
+ f dv 



(« - 2m; - i)G A ?(a; tl )G A «(a; 1 ,)7A + (« - 2uv - i)G A « (x u )G Ae (x v )>yx 



+ 1(1 - 2« - 2 v )G A c(^„)G A «(^ t ,)7A75 - ^G A? (x tl )G A «(x„) 7A 7 5 ] } + 0(twist 3) (59) 

where A. = x' — y' and e = 1 — Hereafter, we use the notations x M = ux' + uy' and Ga m = G Qp A a for brevity. 

We need to multiply this by similar expansion for the antiquark propagator. The product of the two quark 
propagators has the form: 

Trt a 7 a (x'\^\y')t b ^(y'\^\x') = + + (60) 



where 

= * m ~ £) W,yT b ^\(p 2 9ap- PaP 0)-^y e \y') 

fl(2-e,l-e) r(3-2e) Z" 1 



8tt 4 (-A 2 ) 



3-2e 



^[s',s u ](2mA a G / 3A(s«) - 2iuApG a &{x u ) + iA 2 G a p(x u ))[x u , y'\ 



= * B(2 ^! £) ( F A2)3-t / du uu(2A a Ap - A 2 g a/3 ){[x',x u }D x G XA )(x u )[x u ,y'])' 



(61) 



(~A 2 )s 

+ B ^i 6 ^~^ T (tf)2-ll f du([x',x u }[{[-i(^-uu)D x G Xa (x u )A - l -uu(l-2u)D a D x G XA {x u )A p + a^ [3] 
+ i(- -uu)g Ql3 D x GxA{xu)+ -uu{l-2u)g af ]{A- D)D x GxA{x u )}[x u ,y']) (62) 



(f a )^ =- (2A a A f3 -A*g a0 ) B{2 j£ e l J^J^- j\ u dv Tr{uv6(u - v)t b [y',x']t a W,^} 
x GAt(x u )[xu,x v ]G£(x v )[x v ,y'} + uv9(v - u)t a [x ', y']i b [y', £ m ]Ga£ (£«)[£«, x„]G / |(x„)[^,a: / ]| 



B(2-e,-e) r(2 - e) 
I67? (-A 2 ) 2 - e 
1 



Tr{ - 2{A a 5 x + A P 8 X - g aP A x ) (o{u - v)t b [y', x']t a ['x, <B* 



x {(u - 2uv ~ -)Ga£{xu)[x u ,Xv\g£{x v ) + (v - 2uv - -)G^(x u )[x u ,x v ]GAc( x v)}[xv,y'] 
+ 9(v - u)t a [x' ,y']t b [y' , x u ]{(u - 2vu - -)Ga£ x v ]G x i (x v ) + (v - 2vu - ^G^x^Xu.x^Ga^Xv^Xv.x']^ 

- B ^ ~ 6 £ ^ ~ - ^ 2 ~ 2 2 £ ] J dudv Tr { A 2 t a ( [x' , x u ]G ai (x u ) [x u , y'} t b [y' , x v ] G/ (x v ) [x v , x'} + a <-> f3) + 2 [uv + vu) 

x t a {^g a p[x l ,x u ]GAi{x u )[xu,y\t b [y'-,x v \Gl(x v )[x v ,x\ ~ [x' , x u ](GA a (x u )[x u , y]t b [y', x v ]G A pX(v) + a <-> /3)[x v ,x']j 

- 2uA a t a [x', x u ]Gas(xu)[x u , y']t b [y', x v ]Gfi(x v )[x v ,x'] - 2vA a t a [x', x u }G^(x u )[x u , y']t b [y', x v ]G£(v)[x v ,x'] 

- 2uA p t a [x' ,x u ]G A i;(xu)[xu,y'}t b [y' ,x v ]G^(x v )[x v ,x'] - 2vA p t a [x' ,x u ]G a ^x u )[x ul y']t b [y' lXv ]G [{x v )[x v ,x']^ (63) 
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(a) 





(b) 



(c) 




(d) 




(e) 




(f) 





(g) (h) (i) 

FIG. 11: Quark-loop contribution to the gluon propagator in an external field 




(j) 



where we have omitted terms ~ e ct ^ M „A A '(...) l/ and [G A £, Gj] which do not contribute to Eq. looli with our accuracy. 

Next we need to substitute the product ilBT|) into the expression iJSSJ. Since we will integrate the expression 
over x* and (to get U X U^) we can neglect the terms proportional to P a (...)p and {...) a Pp. Indeed, using the identity 



Pa 



1 



P 2 <7 a/3 + 2iG af3 P 2 p P 2 



P 2 .g Q(3 + 2iG a 







we get 



1 1 1 1 

P " P 2 g a . + 2iG a . + O a . ~ P 2 P " P 2 «?. ct + 2iG. Q + a Q PQ " 'P 2 



(64) 



(65) 



As explained in Ref. one can drop the terms proportional to P, since they lead to the terms proportional to the 
integral of total derivative, namely 



[22 2 
/ dx»[oopi, -x»pi]t°[-x»pi, -cx3pi](D.$(-x»pi, ...)) 6 
J s s s 

d 2 2 2 

dx*- — {[oopi, -x*pi]t a [-x*pi, -oopi]($(-x»pi, ...))ab} = 

(XX sis 5 5 5 



Using this property one can rewrite Eq. H67JI in the form 

|i|y')*W|i 



(66) 



(67) 



where 



r 2 (2- £ ) 

47r 2 r(4 - 2e) 



(x|P 2 



T(e) 

-P 2 ) 



' ly)ab ~ Sr^-2e) {A ^$)^ y) l Q du (^^^ G ^)[xu,yT b 



+ ^ r(3 r( 4 e)r( 2 1 £ ) ~ {xl J^r lv) l! dU W'^WKGAfiixJ - uS$G Aa (xu)}[x u ,y'}r b 
(T 2 )% = g aP {\ - 2e) = ^ (x| |y) jf du ([x', x„p^ A (x„)[x u , y']) ab + (T 2 )& 



- 2 r 2 (2-e) , r(l + e 



CTO = *Sa/j(l-2e) 



(x| 



27r 2 r(4-2e) v '(-p 2 ) 1+e 



(68) 



|y) / du/ dv Mv([x',x u ]G A e(a;u)[x„,x„]G A s (x„)[x t ,,y , ]) a '' + (T3)^ 



and O- means "equal up to the contributions ~ P a (...)p and (...) Q P,g". 

Next we expand the propagator (|55|l near the light cone. The first contribution comes from the T\ term which 
represents diagrams in Fig. E] a- g- The calculation yields 
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Ux'dy 1 (x\ 



1 



\ x >r a T?Uy'\ 



i 



■\y) 



bn 



P 2 g. ll + 2iG.^ + O mit ' ' Va ^ u ' P 2 g v . + 2£G„. + O v . 

1 (d\U ab T(e) | 1 M *5 n 1 fl^G?. 6 r(e) ^ 1 



a ' (— p 2 ) e 



}£lv)-^*(2-e>2-«)(*l£{ : 



a ' (— p 2 ) e J p 



}> 



3 2 r 2 (i- e ) 



1 if a 

327r 2 r(4- 2e)7 ~oJ 



ab 



4i 



i[(2 - e)(z - /): - (x - z)l - (z - z')l] + - e [(2- e)(z - y)% - (z' - y)% - (z - z')\ 



+ 2(z - z')l - 2(2 - e) {Z Z '\*_ e + (1 - e) {z - 2 (2 - e) /' ^ + (1 - e)^— ^| 

(ir-z')* (x-z')i (^-2/)* 



2(1-6) 
3-e 
1 



+ 



(z - z') 2 [a; 2+£ + (z- z'): 2 + £ - (* - ^): 2+e - - y); 2+£ 

(x - z)\^ + (z- z')l +e -{x- z')l +t + (l + e) (X ~ Z) * (Z Z ' ] 



( z - y )l^ + (l + e) 



(z - z%(z' /(3 - e)g lj 



(x - z%- e 

• aS r A ? A A 



+ W - y)l +e + (*- *0 



/\l+e 



(z-Oi - * 

+ ^(1 _ e )(3 - 2e )[2(z - z') £ * + (* - «)5 + (* - y) £ * - 2(x - z')l - 2{z' - y)l}\ ( 



lots 
4ttA, 



l-e A | 

e l ^ QS (69) 



In our "external" field the characteristic distances are of the order of width of the shock wave: z*,z^ ~ 

e m y s/m 2 . As we shall see below, the characteristic distances x» and ?/* are ~ e 7 * 2 a/ s/m 2 so we can neglect z* and 
z* in comparison to x* and/or ?/*. The formula (|69[) simplifies to 



dx'dy' (x\ 



1 



1 



1 

8^ 
4? 



P 2 g, il + 2iG^ + 0. M 1 ; iQ/JVS 1 P 2 5t/ . + 2iG„. + O v . 

l_cd\U^_ T(e) -I 1 , iff 3 B(2-e,2-e 
j3 2 I a ' (— p 2 ) e J p 2 



bn 



B(2-e,2-e)(x\\{- 



16tt 2 



(fa 
a 3 



(S) e ^* + ^*^ ("iSt)' '^^""J d \ z * ([ x *> z *]xD i G i .(z,,x ± )[z.,zl] x [z,,y,] x ) ab 



5 2 r(l-e)r(2-e) f^fa /•»• 2 f 2 ' 2 , 



167r 2 r(4 - 2e) 7 a 



-^r / d-z* / d-z^ ([x*,z*] :c G'. I (z*,a;_L)[z st ,z^ a: G. l (z^a; ± )[z^?/» : ] :E ) 



« 6 



1 - e 



X 



1 

3 -e' 



l-e . A | 



The term ~ T 2 coming from the diagram in Fig. Illt i has the form 



/dx'dy' { x \\\ x ')Tf„{y'\\\y) 
J V P 

i f°°S a ,iy, iasy-' ^ r(2- e )r(l- £ ) [ (4-6)(l- £ ) 
^ 2 i ^ WV 4ttAJ r(3-2e) \ e(2 - e)(3 - e) [ * + 1 V) * 



(70) 



(71) 



32tt 2 

2(2-2e + e 2 ) rl-e iasA 2 j_ 



e(2-e)(l-e 2 ) 



4A 2 



=i +e + (-2/)i +e - Ai +£ 



(1 + e)x*y* 



1 



i(2-e)as a2 « 2 5 2 a4 
2A2 Ax ~ 16A! A± 



Ai" e J e(l-e)(2 + e) 
y* s 



ex*y*A 
(2-0(l-e) 



€-2l 



A 2 



[^ 2+£ + (-2/) 2+e -A 2+e -(2 + eKy*A £ J I / d-^([a„z.] ss D x GfA.(«.,a±)[«.,tf.]S 6 



where we have neglected z* in comparison to a;*, as discussed above. Since there is no field outside the shock wave, 
[25| at x*y* > the contribution (|71fl vanishes and at > 0, < one can extend the limits of integration over z* 



18 



to ±00 and get 

jdx'dy 1 { x \±\ x i)T£.{y'\^\y) 



(72) 



x»>0>y„ 



' da I i \ e 



327r 2 J a 3 \asJ V At: A 



l -\i^as r(2- £ )r(l-e) |- (4-e)(l-6) 



er(3-2e) L (2 - e)(3 - e) 



-y)l 



Al — ex*y*A 



1 



,_ 2 2(2-2e + e 2 



1 — e iasA 2 , 



(l-e)(2 + e) 



J (2-e)(l-e 2 ) 
(2-e)(l-e) i(2-e)as 



4A? 



K +£ + (-») 



1+6 _ 



(1 + e)x*y* 



A! 



2A2 



A] 



IQAi 



A 



[** +£ + (~y)l +e - A 2+e - (2 + e)x*y*Ai\} [DG]f 



where [DG] is defined by Eq. I|21|l . Similarly, the integral f^'d^z^dx*, z r ) x D x G\,(z*, x±) in the second term in the 
r.h.s. of Eq. lf7U|) can be reduced to [DG] X . The case a:* < 0,y* > is obtained from (|72*|) by the substitution x <-> y. 
The contribution of diagrams in Fig. Illl i.i is 



jdx'dy' { x \\\x')T£.{y'\\\y) 
J V P 

g 2 , . f(fa/4i\7 

= -|_B(l-e,2- e / — — ( 
lD7r z j Vas/ V 



(73) 



ias 



I (4-e)(l-e) 
Je(2-e)(3-e) 

1 — e iasA 2 ^ 



k + (- y ): - a: - ( Z - zo: 



(4-e)(l-e)x^ 



(2-e)(3-e) (2-e)(3-e)A 



2-e 



A, 4A 2 



2(2-2e + e 2 ) 
<l-e 2 )(2-e) 



*i +£ + (-y) 



l+e _ Al+e 



(1 + e)a:,?/*A* 



(2-e)(l-e) i(2-e)as 



A 2 



A 2 , - ———A 4 



[z 2+£ + (-y) 2+£ - A 2+£ - (2 + e)x*y*A% 



9ij- 



2A2 ~ x 16A4 - 1 

2 - 2e + e 2 



.2-6 



|~1 — e iasA 2 ^ 

at + ^at 



e (l- e 2 )(2- e ) 



e (l-e)(2 + e) J 

[^ +£ + Hy)i +£ - Ai +£ - (i + Oirf 1 ] 



(2-e)(l-e) i(2-e)as 



>,2 ,,2 



A 2 

1 — e iasA 2 
1 i 

A, 



2A2 



4A 2 



a 16At ± 
~2/)i +£ -Ai+ e 



r 2+e 



+ (-y) 2 + £ -A 2 + e -(2 + eKy,A e J 

e(l-e)(2 + e) 9 ^ 

xl +e + (-y*) 1+e - K +e - (i + Oi.j.aj- 1 



2 A, A2 AlAj 



(l-e 2 )(2-e) 



x [x*, z*] am 2Tr{t m G. l (z* 1 ii)[z„ x±)[z^ z*]x + m <-> n}[z»,y»]3 



The final expression for the light-cone expansion of the quark-loop contribution to gluon propagator in the sum of 
the expressions (|70ll . (|71|l . and (|73|l . A very important observation is that the contributions proportional to 

g A n f Jdz*J dzi 9(z - z') (z - z'yG mi (z,)G mi {z',) (74) 

present in the Eqs. I|70l) and l|73[) cancel in their sum. If it were not true, there would be an addtional contribution 
to the gluon propagator Q at the g A level coming from the small-size (large-momenta) quark loop. Indeed, the 
calculations of Feynman diagrams with the propagators JjJ and ljl2|l implies that we first take limit z*, z'^ — > and 
limit dj_ — > 2 afterwards. With such order of limits, the contribution l|74|) vanishes. However, the proper order of 
these limits is to take at first d± — > 2 (which will give finite expressions after adding the counterterms) and then try 
to impose the condition that the external field is very narrow by taking the limit z*,z^ — > 0. In this case, Eq. i|74|) 
reduces to g 4 nf[GG]. The non-commutativity of these limits would mean that the contribution p-[GG]4j- should be 
added to the gluon propagator (JJJ to restore the correct result. Fortunately, the terms ~ (I74|l cancel which means 
that there are no additional contributions to the gluon propagator coming from the quark loop inside the shock wave 
(= quark loop with large momenta). 
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Since there is no external field outside the shock wave, after cancellation of the terms ~ (z — z') e we see that at 
> the sum of Eq. H7U|) and Eq. (|73[1 vanishes, and at > > y„ one can extend the limits of integration in 
the gluon operators to ±oo and obtain 



X* r-X* 



dzj dz:{[x* 1 z*] x G. l {z^x±)[z^z:] x GM*^±-)K,y*\x) ab - [GG]f (75) 

Jy, 

*dz* Hdzl [^z^Tr^G.^,:^^ -> Tr^l^C/]} 



We get 



dx'dy' (x\ 



\„l\ma(rpab , rpab .rpab 

Pig. a + 2iG. a + 0. a lX } ViaP + 'w + ltopMl 



P 2 gp. + 2iG p . + Op. 



\y) 



bn 



x,>o> y , 1 D ,„ „ 1 f djg7 Qb r(e) | 1 

a ' (— p 2 ) £ I p ; 



8tt 2 
zas 



^(2-e,2- e )(a:|-^{- 



32tt 2 



B(l-e,2-e) 



1 (fa / 4i \ e 
a 3 V as 



4ttA 
1 — e iasA 2 i 



1-6 A | 



4A 2 



U e(3 - 2e) 
2(2-2e + e 2 ) 



e(2-e)(3-e)' 

[*i +£ + (-y)i +e - Ai +e - (i + eKi/.a;- 1 ] 



(2-e)(l-e) i(2-e)as 



A 2 

A2-e 



2A2 



A i - 



-A 



[z 2+e + (-y) 2+£ - A 2 + £ - (2 + eJs.y.A;] 



1 — e iasA^ 



A, 4A 2 



16AJ 
2 - 2e + e 2 

Ki-^ 2 )(2- £ ) 



e (l-e)(2 + e) 
H +e + (-y)l +E - Ai+ £ - (1 + e)x*y*Ai- 1 } 



(2-e)(l-e) , z(2-e)a SA2 



A 2 



2A5 



16A1 



A 



„2+e 



(-y) 2 + £ - A 2 + £ - (2 + e)x*y*Al 



e(l-e)(2 + e) 



1 - e 
A* 



iasA^ 
4A 2 



ci +e + (-y)i +e - A^+ e 2 + - A^+ £ - (1 + e^AJ- 1 



e(l + e) A, (l-e 2 )(2-e) 

x\ +t + {-y*) 1+t - Ai+ e - (1 + e^AS" 1 



4Tr{t a d % U x t b d l U x } 



A 2 



AW 



(l- e 2 )(2-e) 



(76) 



We see that the light-cone expansion of gluon propagator contains only Wilson lines and their derivatives as should 
be expected after cancellation of the "contaminating" terms (I74|) . 

Next, to get the expansion of [co,0] s ® [0, — oo] y near the light cone we integrate the expression l|76[l over from 
to oo and over from — oo to 0. It is easy to demonstrate that 



dx, 



dy* ( 



4vrA, 



l-e A | 



(4-e)(l-e) 



1 — e iasA, 



2(2-2e + e 2 ) 
£ (2-e)(l-6 2 ) 



i(2 - e)as_ a2 _ a 2 s 2 ^ 



2A2 



16A4 



4A 2 



^ e (2-e)(3-e) 



< + (-2/): - A e , - ex,y,Al' 2 ] 

1 



e (l-e)(2 + e) 



(2-e)(l-e) 
A 2 



-y) 2+e -A 2+e -(2 + 6)^y,A £ J 



(in particular, it means that the term 172|) coming from the diagram in Fig. Illh does not contribute). The result of 
the integration of Eq. (|76[1 has the form 



[00,0]*® [0,-00]^^^ 



^Ar,B(2-e,2-e)( Xx \^ ± U-^}^\y ± ) + ^rjB(l-e,l 
({ 



r(-l-2e) 



(4-e)(l-e) , 7 + e 



6e(l + e) 6(1 + e) 3(l + e)(2 + e) 



{piYhV 8tt 2 ' v ' (A^)- 1 - 26 



(77) 
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and therefore 



df[oo,0] x ®df[0,-oo] 



x± ^ y± a 2 s ^ B{2-e) r(-2e) 2 b a 2 Ar, T(-2e) 



7T^ 



e (A^)- 2£ ^ x 4tt 2 v y (A 
X ({ - I + 2 - ^} 2^^^} _ 2Tr{t° di U x t%Ul} 

so we obtain 



(78) 



Tr{d?U x dVU}} 
2 



+ 3(l + e)J 



x± ^ y± ajiif B(2-e,2- e) T(-2e) Jja , 
7r 2 e 

La . 4e A » A 

3(l + e)I^ + 3(1 + 6) Ai 



-—rf<>lu?+ £|^A^(i- e ,i- e ) r( " 2,i 



(Ai)- 2 



(Ai) 



2 \-2e 



(79) 



Last, we need to write down the sum of 1/e counterterms to diagrams in Fig. 1)11(1 a-g. It can be read from the 
first term in the r.h.s. of Eq. (|70|l : 



12n 2 e 1 (Ai)~ 2£ x 
Adding the counterterm JHOJ to Eq. lj%T)l we get 



(80) 



TrKt/^t/t}*^ ^Ar? 



B(2-e,2-e) T(-2e) 1 T(-2e) 



(Ai)- 2£ 12e(A 



2 \-2e 



a 2 n f . . , i(— 2ej r/ ^ i \ - 



r(-2c) 



47T 2 iV c 

which coincides with Eq. (|22ll . 



(A 



4e AjAj 



3(1+6), 



3(1 + e) Ai 



Tr{cW^[/t} (81) 
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